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ABSTRACT 

A digital computer simulation program using Digital Simulation 
Language is produced to study the performance of a hydrofoil in calm 
water. Various automatic control systems are studied with the model 
constrained to the pitch-heave-surge mode of operation. 
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TABLE OF SYMBOLS 



a f 


feet 2 


Area of foil 


A s 


feet 2 


Area of strut 


CG 


— 


Center of gravity 


C DF 


dimensionless 


Drag coefficient of a foil 


C DS 


dimensionless 


Drag coefficient of a strut 


C L 


dimensionless 


Lift coefficient 


c s 


dimensionless 


Side force coefficient 


F DiF 


pounds 


Drag force on a foil (water axes) 


F DiS 


pounds 


Drag force on a strut (water axes) 


F LiF 


pounds 


Lift force on a foil (water axes) 


F SiS 


pounds 


Side force on a strut (water axes) 


F x 


pounds 


Force in direction of body X-axis 


F 

XiF 


pounds 


Force on a foil in body X-direction 


F XiS 


pounds 


Force on a strut in body X-direction 


F v 


pounds 


Force in body Y-direction 


F 

YiF 


pounds 


Force on a foil in body Y-direction 


F YiS 


pounds 


Force on a strut in body Y-direction 


F z 


pounds 


Force in direction of body Z-direction 


F 

ZiF 


pounds 


Force on a foil in body Z-direction 


F ZiS 


pounds 


Force on a strut in body Z-direction 


g x 


feet/second 2 


Component of gravity in body 
X-direction 
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9 y 


2 

Feet/second 


Component of gravity in body 
Y-direetion 


g z 


2 

feet/second 


Component of gravity in body 
Z-direction 


h s 


feet 


Height of height sensor above 
instantaneous water surface 


*xx 


2 

Slug- foot 


Mement of inertia about body X-axis 


1 YY 


slug-foot^ 


Moment of inertia about body Y-axis 


r zz 


2 

Slug-foot 


Moment of inertia about body Z-axis 


L 


foot-pounds 


Moment acting to produce roll about 
body X-axis 


^XhS 


& feet 


Distance from CG along X-axis to 
height sensor 


hiiF 


feet 


Distance from CG along X-axis to a 
point of foil force 


^iF 


feet 


Distance from CG along body Y-axis 
to a point of application of foil force 


L ZiF 


feet 


Distance from CG along body Z-axis 
to a point of application of foil force 


^ZiS 


feet 


Distance from CG along body Z-axis 
to a point of application of strut force 


M .. 


foot-pounds 


Moment acting to produce pitch about 
body Y-axis 


m 


slugs 


Mass of the craft 


N 


foot-pounds 


Moment acting to produce yaw about 
body Z-axis 


P 


radians/ second 


Roll rate about body X-axis 


Q 


radians/second 


Pitch rate about body Y-axis 


R 


radians/second 


Yaw rate about body Z-axis 
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S EiF 


feet 


Submergence of a foil in earth axes 


S EiS 


feet 


Submergence in earth axes of the 
assumed point of application of strut 
force 


SHP 


— 


Shaft horsepower 


T 

X 


pounds 


Magnitude of thrust in direction of 
body X-axis 


U 


feet/second 


Velocity in direction of body X-axis 


U E 


feet/second 


Velocity along earth X-axis 


U. 

1 


feet/ second 


Foil velocity in direction of body X-axis 




feet/ second 


Total relative velocity 


V 


feet/ second 


Velocity in direction of body Y-axis 


V E 


feet/ second 


Velocity along earth Y-axis 


V. 

1 


feet/ second 


Foil velocity in direction of body 
Y-axis 


w 


feet/ second 


Velocity in direction of body Z-axis 


W E 


feet/ second 


Velocity along earth Z-axis 


w. 

X 


feet/second 


Foil velocity in direction of body 
Z-axis 


*e 


feet 


Distance along earth X-axis from 
reference origin 


y e 


feet 


Distance along earth Y-axis from 
reference origin 




radians 


Angle of attack of a particular foil, 
designated by subscript 


$1 


radians 


Angle of side slip of particular strut, 
designated by the subscript 


7 ) 


dimensionless 


Propulsion efficiency 
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p 

0 

t 

© 



3 

slugs/foot 

radians 

radians 

radians 



Water density 
Roll angle 
Yaw angle 
Pitch angle 



Subscripts 



C 

P 

s 

M 

F 

S 

i 



Refers to Center foil or strut 
Refers to Port foil or strut 
Refers to Starboard foil or strut 
Refers to Mid foil section 
Refers to foil 
Refers to strut 

When "i" appears in the subscript to a symbol, it 
indicates that the symbol is to be repeated with i 
successfully replaced by C, P, S, and M, as indicated 
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I. INTRODUCTION 



A . COORDINATE SYSTEMS 

Three coordinate systems will be used in this study of hydrofoil 
craft dynamics . They are: the body axes , the water axes and the earth 
axes. Each is a right hand orthogonal system. 

The body axes coordinate system (X, Y, Z) shown in figure 1-1 
has its origin at the center of gravity of the craft and is fixed relative to 

the craft. This system is initially coincident with the origin of the earth 

! * 

axes system, however, it moves with the craft as time progresses. X is 




Orientation of the Body Axes with Respect to the ' 

Craft and Directions of Positive Velocities 
Figure 1-1 

positive forward, Y is positive to starboard and Z is positive 
downward . 

The earth axes coordinate system (X p , Y , Z ) is fixed relative to 

ill Xj £j 

the earths surface. The origin may be chosen at any convenient location, 
provided the body system origin is initially at the same point. The X £ 



11 



axis lies in the horizontal plane; it is initially coincident with, and 
positive in the same direction as, the body X axis. The axis lies in 
the horizontal plane and is positive to starboard when the observer is 
facing. the positive direction. The Z^, axis is normal to the horizontal 

and is positive downward. Water surface motion and crafts' motion 
relative to calm water are described in this coordinate system. 

The water axes coordinate system is aligned with the relative 
velocity vector and resolves into 'lift', 'drag' , and 'side force' directions. 
All hydrodynamic model test data is taken and presented as plots in terms 





Definition of Positive Directions for Euler Angles' 

Figure 1-2 
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Detailed derivations of the transformation matracies between coordi- 



nate systems are outlined in Ref. 1 , however, only the results will be 
considered here. In general, the body axes are displaced from the earth 
axes by the Euler angles, df> , & , and Y • Figure 1-2 defines the posi- 
tive directions for the Euler angles and established the means for trans- 
formations between earth and body axes . 

The hydrodynamic force and moment data are obtained in a water 
axis system which is always oriented with respect to the relative water 
velocity. These data must be transformed into body axes so that their 
effects can be included in the equations of craft motions. In general, the 
water axis orientation is always changing with respect to the body axes. 
This necessitates a transformation to resolve lift, drag, and side force 
quantities into body axes . The orthogonal water axis coordinate system is 
defined in figure 1-3. The angles of rotation between water and body axes 
are c< and and will be defined later. 




Figure 1-3 Definition of Water Axes 
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B. EQUATIONS OF MOTION 



The development of the six equations of motion for a hydrofoil 
follow the same classical lines as those for a displacement type hull. 
This development is clearly outlined in Ref 1 and 4 and therefore 
will not be done here. The six equations for the hydrofoil are: 



Translation 

II = W F X -QW + RV (1-1) 



V= niFylPW * RU 



W~ m f Q U. ** P V 



(1-3) 



Rotation 



P =j xx [l - QR(I„ - I n )+( r +QP)lJ (1-4) 



Q^[n-RP(j; r i 22 ) - (p-R l )i x£ ] 

R = -PQ (lyy-lj- (Q R-P) I xe ] 



(1-5) 



(1-6) 



Comparison of actual and simulated data show that the translation 
terms RV, PV, and PW may be neglected. Similarly, the product of 
Inertia quantity, L^, is less than 10% of 1^ and 4% less than I . 
Consequently, when considered in conjunction with the magnitudes of 
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2 2 

(R+PQ) , (P -R ) and (QR - P) , the terms containing I^ z may be neglected 
for this craft . 

Given the equations 1-1 through 1-6 for linear and rotational 
accelerations in body axes , the velocity terms in body axes are obtained 
by merely integrating each equation one time. These outputs must then 
be transformed to earth axes to be of useful form. The transformation of 
the linear equations become: 

U E = UC OS e Co$Y f V(CosTSin©Sin0-Sinj f Sim (p) 

•fW(cos Y Cos & Cos <p -f- Sin Y Sinj (1-7) 

V £ -/^SiNrCose+V(CoSr Cqs0+ SinYSiWS S/N0) 

d-8) 

+ W(siw YSiho Co s 0 - Cos Y Sim 0) 



— £|Sl/\!©-f'VC-O602if'40-{~ VJ CoS 0 CoS 0 (1—9) 



Only the crafts velocity in the downward direction is required to 
compute the crafts height above the water surface. Therefore, the verti- 
cal position in earth axes is given by: 



Z £ = J"w E d-t 



( 1 - 10 ) 



The crafts position in earth axes coordinates is found from: 

^E“ X Ug dt dt (i-ii) 

The Euler angles ,<0 ,& , and Y , are also of interest and can be 

. 

found by integration of the Euler rates 0 , © , and ^ . The rates in 
turn must be derived from the body angular rates. The Euler rates are not 
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easily obtained because they occur about axes which are not orthogonal. 
This fact can be appreciated by referring to figure 1-2 and recalling how 



the angles were defined. Y is a rotation about the Z q axis, Q is a 
rotation about the Y axis, and 0 is a rotation about the X axis. The 



results in the following set of equations: 



0^ PH- YSim © 

0 = Q CoS 0 -RSiN0 

T “ (q$<n 0H*RCos^Co&©H*^0-*P)SjK © 



( 1 - 12 ) 



d-13) 



(1-14) 



Studies of this craft show that even under the most drastic abnormal 
situations, pitch angles are not expected to reach 10° and under normal 
operating conditions, roll angles will not exceed 3°. Therefore, the 
following assumptions are justified: 



The Euler angles are obtained by integration of equations 1-12, 

13, 14. 

c. calculation of velocity COMPONENTS 

The various foil velocity components must be determined at this . 
point in order to calculate the angles of attack and side slip. The total 
velocity components of each foil in body axes can be expressed in terms 
of the craft linear and angular velocities by considering the craft geometry 
as represented in figure 1-4 and Table 1. 



Sin 0 = Q 
Cos Q = 1 



(1-15) 
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The resultant foil relative velocity equations become: 
Center foil 



U-c. LA + 1— 2c p Q 


d-16) 


Vc-V* - L ic p P+ L xcr R 


(1-17) 


w c = w- L xcf Q 


(1-18) 


Port foil 


C(p ~ U +" Lfpp i-ypp R 


(1-19) 


Vp-V~ L^ppP + L xPF R 


U- 2 °) 


W p = W~LyppP~ L x pp Q 


(1-21) 


Star! oard foil 


U s - U+ - L y S p R 


(1-22) 


= Lz $p P + ^~XSF ^ 


(1-23) 


W s = W+ LySp P “ L xsf Q 


(1-24) 
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Key: 

© = point of application of foil forces 
o = point of application of strut forces 



Note: 

ScCF = ScCS 

Scpf = Scps = SesF = Sess = Scmf 

SfSF = S^SS = _L YPF = _L YPS 
L ZPF = L ZSF = SlMF 

^ps and ^SS are variables and in general will not be equal. 
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TABLE 
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Mid foil 





(1-25) 


Vm “ P *L KMfr R 


(1-26) 


- W- L khp q 


(1-27) 



These foil velocities are in body axis terms and can be used directly 
to calculate the angles of attack and angles of side slip. 

The angle of attack of the foil is shown in figure l-5a and defined 
as the angle whose tangent is the total relative velocity between the water 
and the foil in the body Z direction, divided by the total relative velocity 
between the water and foil in the body X direction. 

<K= ARCTAN -y*: d-28) 

If the foil is not aligned with the body X axis, then the foil has a fixed 
angle of attack and the total angle of attack becomes: 



= ^fUei +c< 



(1-29) 



The side slip angle is shown in figure l-5b and can be represented 



by: 



•£>;. = ARcsin 



V, 



YU 



yiCi + v, i ; + w,\. 



(1-30) 



During actual craft operations , the angles of and will be 
less than 10° permitting use of the small angle approximation. Also, the 
V and W velocities can be considered very small in relation to the U 
velocity.. Applying these simplifications, the equations foro< and 
become: 
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Body X-Axis 




ANGLE OF ATTACK 



Body' X-Axis 




(b) 

ANGLE OF SIDE SLIP 



Figure 1-5 
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(1-31) 



W<L . 

°<i~ U + ^L-VxeJ 

§1 - -¥cL (1-32) 

V U 

D. EXPANSION OF FORCE AND MOMENT TERMS 

With the aid of figure 1-4 and Table 1, the forces and moments for 
the craft are easily written. It is important to note, that when a numerical 
value is substituted for one of the alphabetical dimensions, the sign must 
be positive or negative depending on the relative position from the CG. 

One o r i the assumptions made upon commencing simulation was that 
the struts contributed zero lift and that the foils contributed zero side force. 
Therefore, all body Y motion originates at the struts and all body Z motion 



originates at the foils . 

The force and moment equations become: 

Rr FvcF + F*PF f Rsf + R<mf + F kss +T x (1-33) 

Fy = F yc S +Fyps + F yss + (i-34) 

^2 = F £CF +F ZPF + F 2SF + F Zf1F (1-35) 

L~ (Fg. L y)pp ^(rjLy)^ *" -(F'y L 2 ) ss ~^Fy L 2 )cs (1-36) 

M = ~(F i L x ) CF -(F i L >i )p F -(FzL^sf- -h (F x L 2 ) cp 

SF+(f r x L i )^ F + (FxL?) cs 

^( F x. L l)pS t ( F X l ~ i )ss + T X L ir (1-37) 

N _ ( F y F ( F y l-x)ps ^ ( F y ^x)ss (1-38) 
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All of these forces have a hydrodynamic origin except for the 
thrust and gravity terms. The thrust is always associated with the body 
axis and requires no further expansion. The gravity terms can be simpli- 
fied using the small angle approximation for the Euler angles. (See 
summary of equations A-31, 32, 33.) 

The general expression for hydrodynamic forces in water axes 
coordinates is 

Fc=C^/RCl L=L,0 ) S (1-39) 

where L, D, S represent lift, drag and side force respectively. The lift, 
drag, and side force coefficients vary as functions of angle of attack, 
angle of side slip, submergence, velocity, flap and rudder positions. 

Ideally, mathematical expressions would have been developed to 
correctly depict the interrelation of all the variables which affect the 
force coefficients. However, hydrofoil technology has not advanced to 
a position which would yield such expressions. Some of the major prob- 
lems encountered in deriving a mathematical expression are: 

1) The occurrence of cavitation and ventilation which are 
completely unpredictable. 

2) Even more unpredictable is the cessation of cavitation and 
ventilation. 

3) The nonlinearity of the hydrodynamic coefficients. 

4) Lack of sufficient test data to completely describe all of 
the above . 
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II. HYDROFOIL SIMULATION 



Given the equations of motion, the next step was to formulate a 
program which would simulate the model. The computer language chosen 
was Digital Simulation Language or DSL. DSL was chosen over FORTRAN 
because of its adaptability to control problems and the ease with which 
the user can insert engineering relationships with a minimum amount of 
programming. 

This computer simulation is designed for the HIGH POINT PC(H)-1 
which has a canard foil-strut configuration. Drawings of the craft and 
the foils are shown in figures 2-1 and 2-2. The after foil is one unit, but 
in the equations, it is considered to be divided into three separate segments, 

i.e. mid panel and port/starboard outboard panels. 

The following assumptions were made prior to beginning simulation: 

1 . Craft equations of motion are valid only for the foilborne 

mode. 

2. Weight of the craft remains constant. 

3. Hydrodynamic coefficients are based on fully wetted surfaces, 
i.e. no cavitation or ventilation. 

4. No constraint was placed on crafts heading. 

5. Foils are considered flat surfaces vice dihedral/anhedral. 

6. The craft has no automatic control system. 

7. The craft may be operated in calm water or in a seaway. 

Before the main program could be assembled, subroutines had to be 
written to find the hydrodynamic coefficients of foil lift and drag, and strut 
drag and side slip. Using the simulation curves in Ref. 3 as a guide, the 
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Drawing of Foil and Strut Configuration 



fez £q>l Tqrward £otc 
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Drawing of Foil and Flap Configuration 




I 
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Figure 2-3. Block Diagram of Equation Relationships 



curves were stored as data points in the SETUP routine. The subroutines 



INTERP and INTRP1 are interpolation routines to obtain the proper coeffi- 
cient from the curves. For INTERP, values of angle of attack and flap 
angle are used to enter the curves and obtain a value of foil drag. Angle 
of side slip and submergence are used to enter INTRP1 to obtain strut drag. 
The curve for the lift coefficient of a foil in the fully wetted region is a 
straight line so this coefficient is found by merely solving the equation 
for the straight line. Curves plotted from the subroutines are shown in 
Appendix B. Listings of the subroutines are shown in the Computer Program 
section. The main program was then assembled using figure 2-3 as a refer- 
ence for data flow. A complete listing of the main program is shown in 
the Computer Program section. 

A. SIMULATION OBJECTIVES 

The object of the simulation was to commence the runs with the 
craft foilborne and in a steady state condition. Once this condition was 
achieved, the step-response of the craft would be studied by applying 
step functions to excite motion along and about each body axis . That is , 
perturbations should be introduced to separately excite: 

1. Motion along the X-axis 

2. Motion about the X-axis 

3 . Motion along the Y-axis 

4 . Motion about the Y-axis 

5. Motion along the Z-axis 

6 . Motion about the Z-axis 
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After the step response was obtained, a sinusoidal sea would be 
inserted to observe the crafts motion in a seaway. 

B. PROBLEMS ENCOUNTERED 

There were three major problems encountered during the simulation 
that are worthy of discussion: (1) the unknown stability of the craft, 

(2} specification of units used in the derivation of the equations, and 
(3) numerical quantities for thrust versus speed. 

1 . Craft Stability 

One of the first assumptions made was that the craft would 
be stable in the foilborne mode with no control system . Starting from this 
assumption, the first consideration became that of balancing all of the 
forces and moments. A speed of 36 knots was chosen to insure that cavi- 
tation did not exist at any of the wetted surfaces and that the control sur- 
faces were zero. The fixed angle of attack was assumed zero. 

It was readily apparent that zero fixed angle of attack was 
an invalid assumption because the angle of attack was to be the control- 
ling parameter in stabilizing the model . Reference 5 listed the fixed 
angles of attack as 0 .332 forward and 0.0558 aft. These values were 
used, but the model became unstable at approximately 4.9 seconds. Only 
slight variations of angle of attack above and below a critical value caused 
the model to either raise out of the water or to sink. This pointed out the 
strong dependence of both lift and drag upon the angle of attack and the 
problem of balancing the force and moment equations . Referring back to 
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equations 1-33 through 1-38, one can appreciate the complexity involved 
in balancing these equations as they are all interrelated. In order to 
simplify the problem at all, only the forces in the body X and Z directions 
and the moment about the Y axis were considered . Another simplification 
was that of considering the entire forward foil as a control surface. The 
Angle of Attack program was then written to solve for the angle of attack 
which minimized the force in the body Z direction and the moment about 
the body Y axis . Once the force and moment were minimized , the program 
then checked the force in the body X direction and adjusted thrust to mini- 
mize it. A listing of this program, which is straight FORTRAN, is shown 
in the Computer Program section . The angles of attack determined from 
this program were then inserted into the main program. Significant results 
were obtained, however, the model remained unstable. This indicated 
that there was no feedback to dampen the error, so the model was function- 
ing as an open loop system. Considering the output quantities to be linear 
accelerations and velocities, angular rates and angular positions, while 
the input quantity is reference height above calm water and X„ is the 
measured quantity which is fed back, the block diagram can be represented 
as shown in figure 2-3. No graphs are shovsn of this mode because of the 
instability of the model. 

Stability of the model without a control system was questioned 
at this point and it was determined that assumption #6, i.e. the craft has 
no automatic control system, was invalid. Reference 5 confirmed the fact 
that the model was only stable in three degrees of freedom without a control 
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Figure 2-4. Block Diagram Of Sysiem With Open Loop 



system. Additionally, it was determined that when the model was con- 
strained to the pitch-heave-surge mode of operation it was still only 
marginally stable without a control system. At this point, a control sys- 
tem was added and it will be discussed in detail in Chapter III. 

2 . Specification of Units 

Considerable confusion resulted from the interpretation of 
the various references concerning the units used in defining Euler angles 
and angular rotations. This problem has been alleviated here by the 
specific definitions in the Table of Symbols. 

3 . Thrust vs Speed 

Since no graphical or numerical data was available to show 
the thrust vs speed relationship, the computer was utilized to determine 
a value of thrust for a given speed. Once the Angle of Attack program had 
found the proper angle of attack to minimize the force and moment, the pro- 
gram then checked the force in the body X direction and adjusted the thrust 
to minimize this force. This produced an initial condition for thrust which 
was then inserted in the main program. The main simulation program incor- 
porated a thrust trim routine which is shown on the main listing. 

The thrust trim routine compares the magnitude of the force 
in the X direction with the thrust and then tests the difference with an 
arbitrarily small fixed error. If the difference is greater than the error, 
then thrust is set equal to the value necessary to make the force in the X 
direction equal to zero. At the same time, the acceleration in the X 
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direction is also set equal to zero. If the difference is less than 



the error, then the thrust is held constant and the acceleration in the X 
direction is allowed to vary. 



% 
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III. AUTOMATIC CONTROL SYSTEM 



Once it became apparent that the model would remain unstable 
without a control system, the only alternative was to add the control sys- 
tem. As was previously stated, the problem will be greatly simplified by 
constraining the model to the pitch-heave-surge mode of operation. In 
so doing, the control system can be limited to the forward foil. It was 
assumed that once the control system for the forward foil was operational, 
it could be easily extended to the after foil. 

A. HEIGHT SENSOR 

For this control system, a height sensor was incorporated which 
sensed the change in height of the CG above calm water. This change 
was utilized as a positional feedback signal to control the flap angle on 
the forward foil. This was accomplished in the simulation by defining the 
reference height of the CG for the particular speed and subtracting the 
actual height to obtain a difference called DELH. This difference was 
multipled by a gain factor to produce the variable FLAP defined as: 

FLAP = DELH*COEFH (3-1) 

Figure 3-1 shows the contents of the automatic control block in the pro- 
gram listing. Figure 3-2 is the block diagram of this system showing 
submergence height as the measured quantity being fed back to the 
controller. 
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The amount of flap is governed by the unbalance of the equations 
of motion. Since the flap angle is directly proportional to both lift and 
drag it makes a significant contribution to the balancing of the equations 
of motion. A parameter study was made and determined that the optimum 
value for COEFH was 0.85. Figures 3-3, 4, 5 show the response of the 
model to a positive 5 deg/sec disturbance applied to Q for 0.2 seconds. 

It must be pointed out that with the configuration used, the flap was 
capable of moving instantaneously because an inertia term had not yet 
been added . 

With the model stable using position feedback, the next refinement 
was that of utilizing acceleration, velocity and position signals for feed- 
back. It was realized very quickly that this was an unrealistic approach. 

B. RATE FEEDBACK 

With the model stable using position feedback, the next refinement 
was that of utilizing acceleration, velocity, and position signals for feed- 
back. It was realized very quickly that this was an unrealistic approach 
in that two of the three feedback signals were measured at the CG, 
therefore, neglecting any pitch component. The equation for flap then 
becomes: 

FLAP = GAIN1 DELH-GAIN3 W-GAIN4 W (3-2) 

where GAIN1, GAIN3, and GAIN4 are gain constants. Figure 3-6 shows 
the contents of the automatic control block for this system and figure 3-7 
is the block diagram. The model became unstable for all values except 
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Figure 3 - 1 . Automatic Control Block for Height 
Sensor 
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Figure 3-2. Block Diagram Showing Height Feedback 
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Figure 3-6 . Automatic Control System for Rate 
Feedback 
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Figure 3- 7 . Block Diagram Showing Rate Feedback 



zero, of GAIN3 and GAIN4. No graphs were plotted because of the 
unstability of the system. 

C. RATE FEEDBACK WITH ACCELEROMETER 

This system was implemented by using an accelerometer to sense 
motion at the forward foil. The accelerometer corrects the body vertical 
acceleration by the pitch rate translated to the forward foil. This pro- 
duces the true acceleration at the forward foil. The equation is given by: 

A ZC - (F Z /m) - I XCF 0 (3 - 3) 

Figure 3-8 shows the contents of the automatic control block for this 
system and figure 3-9 shows the block diagram. This control system pro- 

I 

vided stable operation as long as GAIN3 was less than 0.02 and GAIN4 
was less than 0.2. Figures 3-10, 11, 12 are graphs of. the system 
response to the same disturbance stated previously. 

All three of the systems listed so far were not physically realizable 
in that the flap was provided instantaneous response to any error signal. 
This is not the case in a real system because the flap possesses a finite 
amount of inertia which produces a time lag in the response. 

D . RATE FEEDBACK WITH REAL POLE 

The instantaneous response problem was taken care of in this system 
by adding an inertia term. This entails adding a pole to the root locus 
and placing it on the real axis in the left half plane. The addition of the 
pole helped to stabilize the system, but also increased the overshoot and 
and time lag. The time lag associated with the real pole caused the 
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Figure 3- 8 . Automatic Control Svstem With an 
Accelerometer 
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Figure 3- 9 . Block Diagram Showing Rate Feedback With an 
Accelerometer 
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motion of the flap to lag the motion of the model and thus cause the flap 
to continue moving after the error had ceased. It was also noted that 
when the acceleration and pitch angle gain constants became too large, 
the system went unstable. Comparison of figures 3-13, 14, 15, with 
3-10, 11, 12, confirms these statements. 

The contents of the automatic control block are shown in figure 
3-16. In the block diagram of figure 3-17, the real pole term is incor- 
porated into the controller block. 

E . RATE FEEDBACK WITH COMPENSATOR 

In order to minimize the effect of the inertia term , a compensator 
had to be included in the system. The inertia term caused a time lag in 
the response, therefore, a lead compensator had to be added to make the 
overall response approximate that of the instantaneous response. The 
term that was actually inserted was a computer artificiality because it 
took on the role of both real pole and compensator. The compensator of 
the simulation takes on the form (P ^S+lJ/^S+l) . In the root locus , 
another pole was added, however, it is so far to the left that it was 
ignored . 

Comparison of figures 3-18, 19, 20 with figures 3-10, 11, 12 
shows that the compensator did reduce the overshoot and time lag associ- 
ated with the real pole. Close observation shows that the response of the 
model with the lead compensator closely approaches the response of the 
model with an instantaneous flap. 
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Pitch Angle vs Time 
Real Pole Included 
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Figure 3- 17 . Block Diagram Showing System With Inertia 
Included 
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The contents of the automatic control block are as shown in the 



main simulation program in the computer program section. The block dia- 
gram for the system is shown in figure 3-21. 
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Figure 3-2 1 . Block Diagram for Compensated Control System 



III. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 



An invalid beginning assumption, i.e. the craft is stable without 
a control system, caused the loss of a great deal of time early in the 
research. This was not a complete loss, however, as it necessitated 
the writing of the ANGLE OF ATTACK program to determine the angles of 
attack and thrust which balanced the equations of motion. 

The data obtained from the ANGLE OF ATTACK program illustrated 
the tremendous effect slight changes in angle of attack make on the 
motion of the model. Since these motions are so slight and occur in such 
short time intervals, it would not be possible for manual control. This 
necessitates an automatic control system that will provide the necessary 
outputs to operate the model according to the given inputs . 

A stable model was obtained as a result of this research, however, 
it is somewhat limited. Using the present model and expanding it to 
include the after foils, the model would be completely stable in six degrees 
of freedom . 

Recommendations for further study include; 

1 . Incorporate a complete automatic control system which would 
include all three foils . 

2. Produce a simulation for the transition period between hull- 
borne and foilborne conditions . 

3-. A study of the step response of the model in all modes. 

4. A study of the stability of the model in both regular and random 

seas . 

5. A study of the turning characteristics. 
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APPENDIX A 



SUMMARY OF COMPLETE EQUATIONS 



E quations for Accelerations in Bocp f Axes 

~ rnF x + RV~ QW 

V = m r y t P\V ~ Ru 

V - m F 2 4- Qu “ Pv 

P - 4 [L-Qf«I M -I„)+(*+VP)l xi ] 



( A- 1 ) 
( A-2) 

( A-3) 

CA-4) 
(A-5) 
( A-6) 



Velocities in Bod - ' 7 ' Axes 



U =Jddt 


(A-7) 


V =Jvdt 


(A-8) 


W = 1 wdi 


(A-9) 


j 


P=|Po-i 


( A- 1 0 ) 


Q=/<?di 


( A- 1 1 ) 


R=j r Rdt 


CA-12) 



6i 



T ransformation of Velocities from Body to Earth Axe s 

Uc-U Cos 0 Cos Y +V(Cos Ysm © Sin 0 - Sin Y Coscp) 



Y W (cos yS>in 9 Co$(/> f Sim Y Sin <j>) (A-13) 

Vc = U CsS © SlN^ rV(CoSY COS0 Y SlNY SlNQ S)N 0) 

•bVl(SlN YSlhie CoScp -COSY SlH 0) C A— 14) 

w £ - - usiM© +vcds e sin0 + wcos© c$s 0 ( a-i 5 ) 

$ ~ p + ^ Sin 9 • ( A- 1 6 ) 

©= QCOS0-R$lN0 . (A-17) 

S' - (q S ; k) 4> f RCos 0) CoS & + (0-P) $>in B (a-18) 



Positions in Earth Axes 



X 

m 

n 


Ju E dt 


(A- 19) 




[ Vu-dt 


( A-20) 


II 

Ul 

riJ 


I' W £ dt 


( A- 2 1 ) 


^=j 


f 0 di 


( A-22 ) 


e=J 


' © di 


(A-23) 


y=( 


• 

fdi 


( A- 2 4 ) 
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Kx pansion of Force and M oment Equation s 



F x ^ Z F xiF +S F x i s + mg x -t T x 


(A-25) 


Fy = Z Fy d ' S +• KP^y 


(A-26) 


^ = Z ZZf f • 


( A- 27 ) 


L - (Fgly )p F f (f r iLy) SF “(py/-2)p S _ (FyL 2 ) 5 5"(FyL 2 ') cs 


(A-28) 


/I = Z (f,z. x ) 4> +Z (f x l 2 \y +Z (F x 1 ,Zs + T x L 2T 


(A-29) 


N - Z (FyLjis 


(A-30) 






Expansion of Gravity Terms 









™gx - 


Sin © 


(A-31) 




yng Cos © Sin <p 


(A-32) 


w 3* = 


vr )Cj Cose Cos <j> 


• (A-33) 



Foil Velocity Components 


in 


Body Axes 




Center Foil- 










«c = 




+ L 


2rCF ^ 


( A-34 ) 


V c = 


v- 


^2CF ^ F F X CP F 


( A-35) 


W c - 


w 


“ L, 


:CF Q 


(A-36) 
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Port Foil- 



Up - U. + L iP cQ-L >PF r 


( A-37 ) 


* Vp = V - Li pp P+ L x pp R 


(A-38) 


V/p - W + Lypp P~ Lxpp Q 


(A-39) 


Starboard Foil- 




U S — Ci'h L Q — Ly£ p R 


(A-40) 


Vc,= V- l- 2 5 F P-f-L XS p R 


u-4i) 


— W 4- Ly 5 c P~ L p Q 


(A-42) 


Mid Foil- 




u M =a+i iMF Q 


(A- 43) 


V h - V- L 2A)P P f l xmf r 


( A- 4 4 ) 


W M = Vl'L K XrQ 


(A-45) 


Transformation of Water Particle Orbital Velocity 
from Earth to Body Axes 




~ ^ Cose -Wg W4 ' 


(A-46) 


V W4 = U EW l(ccs tSbV©$/M0 - Sim t Cos 6) + V £W ; (sin ? 


( 


■ S/A/ 8 S/M0 -hCosfCcS <f) + VJ Bw l Cos e 


(A-47) 
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y w i - ^ewI ( as ? $ lN s/w + Vcwi ( SiN ' yS,Nd CcS 4 

- Cos ^S/KJ0) -f- V/ £W [, Cos (j) (A-48) 

Each of these equations must he repeated for each foil. 



Relative Velocity Components 

U r l = ( A-49 ) 

V'Ci - Vi ’ (A- 50) 

Wn=W;-W w ; (A . 51) 

V 

One set of equations for each foil. 

Angles of Attack and Side Slip 

°C L => 1RCTAN W "/' Ul .; + fixed (A- 52)* 




AttCSlM 
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( A- 53 ) * 



Total Relative Velocity at a Particular Foil or Strut 

*2f ~ \/ ul + V* +Wc • ( A— 54) 



* One set of equations for each foil. Angles of attack 

, I 

and side slip are calculated in radians. 
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Hydrodynamic Forces in Wa+er Axe s 

(Each foil equation and strut equation must be repeated 
for each foil and strut, respectively) 



Foil Lift 







(A- 56) 


Foil Drag 






- 


F f>p = H Ap Cdf 


(A- 56) 


Strut Di^ag 




• 




Fqs “ Z A3 Cos 


(A-57) 


Strut Side Force 






F s - i c s 


(A- 58) 


Transformation of Foil Forces from Water Axes 


to Body 


Axes 
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APPENDIX B 



Fwd Foil Lift Coefficient vs Ancle of Attack 




X Scale =4.0 degrees/inch 
Y Scale =0.40 units/inch 
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Mid Foil Lift Coefficient vs tine 




X Scale =4.0 degrees/inch 
Y Scale = 0.40 units/inch 
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Aft Foil Lift Coefficient vs Angle of Attack 




X Scale = 4.0 degrees/inch 

Y Scale = 0.40 units/inch 
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HACT=-Z E+L XH^THE T A 

0 ELH=HR EF-HACT 

AZC = FZ/MASS-LXCF*QD0T 

FL API = ( GA IN1*DFLH+GAIN3*AZC-GAIN4*THETA) ^GAINC 
FLAP = tEDLAG(0 .0,0. 10, 0.16, FLAP 1 } 
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ANGLE OF ATTACK ROUTINE 
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SUBROUTINE SETUP 
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OEFCX( J,K )=INFO( I ) 
COUNT = ICQUNT +1 
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